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Abstract 

We study the spin thermalization, i.e., the inter-spin energy relaxation mediated by electron-electron scattering in small 
spin valves. When one or two of the dimensions of the spin valve spacer are smaller than the thermal coherence length, 
the direct spin energy exchange rate diverges and needs to be regularized by the sample dimensions. Here we consider 
two model systems: a long quasi-lD wire and a thin quasi-2D sheet. 

Keywords: A. magnetic films and multilayers, D. thermoelectric effects, D. heat conduction, D. inter-spin heat 
exchange 
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1. Introduction 

The thermoelectric response of a ferromagnet| normal 
metal|ferromagnet spin valve [l], 0, Q depends sensitively 
on the strength of inter-spin energy relaxation (spin ther- 
malization) inside the normal metal spacer In large 
structures at high temperatures, spin thermalization is 
dominated by electron-phonon coupling, whereas at low 
temperatures direct spin-flip scattering becomes impor- 
tant. A third mechanism is the electron-electron scatter- 
ing, which in relatively large spin valves is weak and can 
typically be neglected. However, for smaller spin valves 
the electron-electron interaction becomes stronger. When 
one or more of the dimensions are smaller than the ther- 
mal coherence length £t = y/hD/ (JsbT), the kernel of the 
electron-electron collision integral should be calculated for 
reduced (ID or 2D) spatial dimensions. In these cases the 
thermalization rate formally diverges 0, Q and needs to be 
properly regularized. In this paper we discuss such regular- 
ization schemes and calculate the resulting thermalization 
rates. 



2. Theory 

In a biased spin valve with an antiparallel configura- 
tion of the magnetic reservoirs, the electron distribution 
function may depend on the spin index a. In this case the 
electron-electron interactions can be described by the three 
different collision integrals represented by the diagrams in 
Fig. [T] and calculated as explained in the Appendix. 
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Figure 1: The three types of electron-electron scattering vertices in a 
system with spin-dependent distribution functions f tT (E). a) Equal- 
spin scattering that does not lead to spin thermalization, b) spin 
conserving scattering and c) spin exchange scattering. The latter 
two cause spin thermalization and are discussed in this paper. 



The kernels of the collision integrals depend on the 
Fermi liquid triplet parameter F > — 1. The precise value 
of F is important for the strength of spin thermalization 
caused by the electron-electron interactions: as shown be- 
low, close to the Stoner instability at F = — 1 spin ther- 
malization is quite strong, whereas it is much weaker for 
FfaO. 

The inter-spin energy exchange can be described by 
calculating the spin thermalization heat current 



Q = VpQ, J dEEI n (E), 



(1) 



where 1'^ is the collision integral for scattering between 
spin f and spin J, electrons, vp is the density of states 
at the Fermi level and f2 is the volume of the spin valve 
spacer. 

In the following we first discuss the resulting inter-spin 
thermalization due to the regular spin conserving term. 
This can be done for an arbitrary spin heat accumulation 
T-f — T|_, where T^-^ is the effective electron temperature 
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for spin f (J,) electrons. Then we turn to the spin exchange 
term, concentrating on the linear response regime in which 
Tf - T|. < (T-f + Tj/2. Within the calculation, we limit 
ourselves to the case in which the difference between the 
chemical potentials of the two spin species, i.e., the spin 
accumulation in the spacer /i s = ^ — [i^ w 0. This is 
because in the limit fi s <^ T the effect of a finite fj, s on the 
inter-spin relaxation is quadratic in /i s . 

For the electron-electron interaction, the effective di- 
mensionality (0D, ID, 2D or 3D) of the spin valve island 
depends on the ratio between the island thickness d, width 
W and length L to the thermal coherence length £t (with- 
out losing generality, we assume d < W < L). Alter- 
natively, we can compare the temperature to the Thou- 
less energies hD/(L, W, d) 2 defined by these size scales. 
Therefore, Od case is realized when fc^T <C hD/L 2 , the 
fD case when hD/L 2 < k B T < hD/W 2 , the 2D case 
when HD/W 2 < k B T < hD/d 2 and the 3D case when 
fc B T> hD/d 2 . 

3. Spin conserving term 

We assume that the spin-dependent electron distribu- 
tion functions can be written in terms of the Fermi-Dirac 
functions with spin-dependent chemical potential fi a and 
temperature T a 0, i.e., f a {E) = f°(E; /v,T ff ) = {1 + 
exp[(E — Ha) / (ksTa)]} -1 . In this case, we can use the 
standard relations^ 



(2a) 



/°(£;/x,T)(l-/°(.E + x; M ,T)) 
= [f(E + x; /i, T) - f{E; fi, T)]n(-x; T), 

J dE[f°(E; /i, T) -f(E + x- /i, T)\ = x, (2b) 



where n(x,T) = {1 — exp[x / (ksT)]}^ 1 is the Bose distri- 
bution function. In the following, we use the short-hand 
notation n„{x) — n(x;T a ). Using these relations, the spin 
conserving collision integral can be written in the form 



coth 



2k B T 



— coth 



2k B T- 



(3) 



Here we used the symmetry of the kernel /Cf,(u;) = JCb(~^), 
by which the terms proportional to the chemical potentials 
/io- vanish. With a kernel K-b(uj) = Kb\uj\ a , this yields for 
the spin conserving inter-spin heat current 



Qb 



VF SlK b T{A + o)Li 4+a (I)4 +a trrA+a 4+ 



x We set !i = ks = 1 in the intermediate results, introducing them 
only in the final results. 



where T(x) is the gamma function and Li n (x) is the poly- 
logarithm function. Using the JCb(uj) from Eqs. (|A.3|i — 
(|A.5|) . we get 
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where (,{x) is the Riemann zeta- function and the super- 
script of Q indicates the dimensionality. This term is reg- 
ular, but can in most cases be neglected since in 3D wires 
the spin thermalization due to electron-electron interac- 
tions is typically weaker than either the direct spin-flip 
scattering (at low temperatures) or electron-phonon scat- 
tering (at high temperatures) [4||. 

4. Spin exchange term 

For the spin exchange term, the collision integral can 
be simplified to 

h=\f dudE% c (u)(fr(E - u) - f t (E' - to)) 



x (h(E') -h(E)) 



. . E' — E\ ( E' — E 

coth — — coth 



2k B T 



2k B T, 



2[n f (E'-E)n^E-E')-7i t (E-E')r H (E'-E)] 

(6) 



leading to the inter-spin heat current 



• v F Vt 
Qc = — ;r-x 



dEEdE'dujK, 



: Hsinh(^Hf^l) 



cosh (|±h) cosh cosh (^) cosh (^) ' 

(7) 

This integral is more difficult to handle than Eq. ^ . There- 
fore, we concentrate on the linear response limit Q c = 
K c (T t - T±) and get 

i^f^^ f dxdydwlC c (Tw)x(x — y) 



K r = ^—T 



16 J cosh (f) cosh (f) cosh (££s>) cosh 

' (8) 

The integral over x and y can be carried out analytically, 
leaving 

v F nT 3 f°° w 2 IC c (Tw)(4n 2 +w 2 ) 
K c = — — — / dw -rr-T- ■ (9) 
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sinh 2 (f ) 



In the 3D case with K. c cx \to\ 1 / 2 , this integral is regular 
and yields the result presented in Ref. However, in 



2 



the ID and 2D cases the integral over the first term has 
an infrared divergence and has to be regularized. In the 
following we present a simple regularization scheme based 
on the finite size of the sample. 

Note that in Eq. (O the cutoff scheme needs to be in- 
voked only for the first term. We include also the second 
term in the 2D case, where the divergence is only logarith- 
mic. 

4-1. ID wire 

As discussed in the Appendix, in finite systems the 
kernels actually are sums over momenta. In the ID wire 
the sum over the quantized momenta q n = 2im/L, where 
L is the length of the wire, can be carried out directly. In 
this case the characteristic scale for the frequency is given 
by the Thouless energy E T D = D/L 2 . When \lo\ > E T D , 
we recover Eq. (|A.3[) . For \ui\ <C EjP, the kernel does not 
depend on oj, 



K. c (\lo 



n 2 F 2 



180n(E^P) 2 {l + F) 



(10) 



We can hence use E T B as a lower cutoff in Eq. In 
principle, we should also include the part of the kernel 
which is constant for low frequencies oj < E^P. However, 
the contribution to the integral from this part is of the 
order of {E T D /T) 3 / 2 times the contribution of the second 
part, and can hence be disregarded in the ID limit E T B <C 
T. The remaining integral is K c oc I\ with 



h 



wdw 



ip it sinh (w/2) 



(11) 



In the limit k B T ^> E T D the most significant contribution 
comes from frequencies w <C 1 for which sinh(w) ~ u>, 
leading to 

/ 

(12) 




Including this contribution, we obtain for the spin ther- 
malization heat conductance from the spin exchange con- 
tribution 



K. 



ID 



F 2 k B ^^2-K(k B Tf 



3HE T D {F- 



2) (f + VfTT+i) 



(13) 



where g = Gq / (e 2 / h) is the dimensionless conductance of 
the contacts. Let us estimate the characteristic temper- 
ature and the thermal coherence length £t- The latter 
is 
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(14) 
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(15) 



Wires with lateral dimensions less than £t fall into the 
ID limit. Connecting such wires to reservoirs via contacts 
with resistance 10 £1 and assuming F = —0.3, we get 



D 



1 /im\ Gq 



0.001 m 2 /s V L J 0.1 S 



(16) 

For wires fabricated from materials close to the Stoncr in- 
stability, such as Pd, say with F = —0.9, this characteristic 
temperature is 
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These values should be compared to the characteristic tem- 
perature due to the electron-phonon interaction using the 
results from Ref. 141. This is 



10 9 Wm~ 3 K- 5 0.015 (/im) 3 G 

0.1 S 



-1/3 



Tch.c-ph ~ 4 KX 

\ Zj II U.±0/ 

(17) 

In this case electron-electron interaction in systems close 
to the Stoner instability is the dominating spin thermal- 
ization mechanism. 

Another way to characterize the spin energy exchange 
is via the spin thermalization time r s t = d^e 2 vpT^lj '(2K C - C ). 
From Eq. ([13]): 

1D _ tt( f + 2) (f + VfTT + 1) hv F n e t u 

Tst '°- C " 8V2F 2 k^f' 

This scattering time should be compared to the direct spin- 
flip time of roughly 100 ps in typical samples 0. With 
F = —0.3 and some other typical values for metals with 
wire cross section A — Wd, r st! c-o is rather long, 



.ID 



F=-0.3 



where we reintroduced k B and h. The diverging term 
hence gives rise to a temperature dependent scaling as 
cxT 2 . 

As discussed in Ref. [4| , a spin valve can be character- 
ized in terms of the temperature above which spin ther- 
malization is stronger than heat diffusion through the con- 
tacts, i.e., K c CqGqT, where Go is the spin-averaged 
contact conductance and Cq = 7r 2 fc 2 3 /(3e 2 ) is the Lorenz 
number. This characteristic temperature for electron-electron 
interaction in a ID sample is 



1 ns 
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1 /xm 



D 



IK 



10 47 J-W 3 1500 nm 2 L 0.001 m 2 /s T 
whereas for F = —0.9 it is 

F=-0.9 „„ 

20 ps 
v F A 



r 1D 

st,c— e 



1 /im 



D 



IK 



10 47 J-^" 3 1500 nm 2 L 0.001 m 2 /s T 

(18) 

These examples show that the electron-electron interaction 
in ID samples is especially relevant for systems close to the 
Stoner instability. 



4-2. 2D square sheet 

In the two-dimensional case K, c 



K2£>M 



and the 



resulting integral has a logarithmic divergence. Now the 
regular term in Eq. @ yields 

and for the diverging term, we need to evaluate the integral 
wdw 



h 



e t /t sinh (w/2) 



= 2^coth(^] -21n 
T \2T 



(19) 
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Here the lower cutoff Et — hD/W 2 is the confinement 
energy due to the smaller of the two extended dimensions. 
Hence, for the full spin thermalization heat conductance 
we get 



2D F 2 k%T 2 ( 47 r 2 (ln(^)+l)+6C(3) 



YlixhE^ (F 2 + 3F + 2) 

where E 2 ? = hD/{LW). 

In this case the characteristic temperature is 



ch,e— e 



(F + l)(F + 2)g 



(20) 



H/rp 



4F 2 W 



(F+l)(F+2)ge~ 



1W 



(21) 

where W{x) is the principal solution for z in x = ze z . Here 
we included also the regular terms from Eqs. (|5bl) and ([2"U1) 
besides the term proportional to ln(T/F/r)- To estimate 
this characteristic temperature, we consider a square bar 
of area (width x length) A — 1 (/im) 2 and much thinner 
than £t- In this case we have T 2 ® c _ c m 10 K for wires 
with F = -0.3 and T 2 h D c _ c w 0.3 K with F = -0.9 and 
otherwise similar values as in the above ID case. For com- 
parison, with the thickness d = 30 nm, we get for the 
phonon contribution T c h. -ph ~ 3.2 K. 
The spin thermalization time is 



2D 



4tt 3 (F + 1)(F + 2)v F nh 



H/rp 



8n 2 F 2 



h ( W 



1 +3(2 + F + 4F 2 )C(3) k b 



(22) 



Assuming a film thickness of 30 nm, for F = —0.3 we get 
Tst,c-o ~ 7 ns whereas for F = —0.9 r stiC _ w 80 ps at 
T = 1 K and otherwise similar conditions as above. 

5. Conclusions 

One of the key features enabling the success of spin- 
tronics in the past two decades is the long spin-flip time 
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Figure 2: (Color online): Temperature dependence of the spin ther- 
malization heat conductance due to electron-electron interaction for 
a wire with length 1 /im, width 100 nm and thickness 10 nm, diffu- 
sion constant D = 0.03 m 2 /s and Fermi liquid parameter F = —0.3 
(solid lines) and F = (dashed lines). The three regimes plotted in 
the figure are for the ID, 2D and 3D limits. The heat conductance 
is written in units of kgE^P /H, where E^P /kg = 22 mK, corre- 
sponding to the crossover temperature between Od and ID behavior. 
In the case F = the spin exchange contribution vanishes, which 
changes the temperature dependences of the heat conductances in 
the ID and 2D cases. 



found in metals. One recent trend has been the study 
of thermal transport in spintronic systems. Special in- 
terest has been devoted to the study of magnetothermal 
effects, where the heat transport depends on the mag- 
netic configuration of the system. Besides the long spin 
relaxation time, such effects rely on relatively weak inter- 
spin energy relaxation. This was studied in detail by the 
present authors in Ref. It was found, among other 
things, that making the sample smaller reduces the effects 
of inter-spin relaxation. However, the conventional 3D the- 
ory of electron-electron interaction in disordered metals 
applies only when the system dimensions exceed the ther- 
mal coherence length Moreover, it was pointed out 
in Refs. 0, Q that the electron-electron contribution to 
the energy relaxation in reduced dimensions diverges. In 
this paper we have addressed this divergence by introduc- 
ing a natural cutoff emerging from the sample size. The 
main results of this consideration, Eqs. (fTB")) and (|2"0)l . show 
that indeed the electron-electron interactions in ID and 
2D samples are more relevant for the spin thermalization 
than in the 3D case. However, for typical normal metal 
spacers for which the Fermi liquid triplet parameter F is 
far from the Stoner instability F = — 1, this effect is still 
masked by direct spin-flip processes at low temperatures 
and electron-phonon relaxation at high temperatures. The 
situation is different for materials close to the Stoner in- 
stability, for example palladium wires doped with nickel. 
In these systems the inter-spin energy relaxation due to 
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electron-electron interaction may be the dominating re- 
laxation mechanism. 

The behavior of the spin thermalization heat conduc- 
tance due to electron-electron scattering is illustrated in 
Fig. [21 It shows the regimes of different dimensionality for 
an example system. 

The theory presented in this manuscript is based on 
the collision integrals calculated in Ref . @ , which general- 
ize the Altshuler-Aronov theory Q to the spin-dependent 
case. We point out that the electron-electron energy re- 
laxation described by this theory in the ID limit has been 
investigated experimentally in detail However, accord- 
ing to these experiments the measured energy relaxation 
is stronger by roughly an order of magnitude than that 
predicted by the theory. This is evidence that the re- 
sults for the spin thermalization conductances presented 
in Eqs. (f5a|) . (|5b|l. (fl"3|) and (|20| might be underestimated. 
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Appendix A. Collision integrals and kernels 

The collision integrals describing electron-electron scat- 
tering in spin-dependent systems are 

I a (E) = J dujdE'JC a {uj)x 

[(1 - U(E))(1 - U{E'))U{E - u)U{E' + W ) (A ' la) 

- U(E)f a (E')(l - U(E - - f a (E' + U ))] 

I b {E) = J dujdE'K. b (uj)x 

[(1 - f a (E))(l - f- a (E'))U(E - co)f^(E' + u>) 

- U(E)f^(E')(l - f a {E - w ))(l - f. a (E' + u))] 

(A.lb) 

I C (E) = [ dujdE')C c (uj)x 



[(1 - U(E))(1 - J_ a (E'))f a (E' - w)f-„{E + w) 
- U(E)f^(E')(l - UE' - W ))(l - f-„(E + u>))]. 

(A.lc) 



The kernels Ka^^u) depend on the dimensionality of the 
sample @: 



ICa = 
JCb = 

K. r 



E 



i 



(I +F ) + (2Dq 2 ) 2 



1-kvfVL ^ uj 2 + (Dq 2 ) 2 (1 + F) 2 



(Dq 2 ) 2 



E 



4 ' (2Dq 2 ) 2 



2Trv F n ^ uj 2 + (Dq 2 ) 2 (1 + F) 2 + 



(A.2a) 
(A.2b) 



(Dq 2 ) 2 



E 



F 2 



2nVF n ^u 2 + (Dq 2 ) 2 (i + F) 2 + (" ( ff a ' 

(A.2c) 

Here D is the diffusion constant in the spacer and F the 
Fermi liquid interaction parameter in the triplet channel. 

In this paper we describe systems with at least one 
extended dimension, i.e., quasi-lD, quasi-2D or 3D metal 
islands. According to the standard prescriptions for cal- 
culating the sums 

/>oo 

V i y L/(tt) / dq, d= 1 

q J ° 

/•CO 

V i— > A/(2tt) / qdq, 

V ^ n/(27T 2 ) / q 2 dq, 
JO 



d = 2 
d = 3, 



we get for d = 1 

4F(F+1) 



K a 
Kb 

K. n = 



(F+2)(\/F+T+l) 



1 



1 



2irAv F ^j2D(F + \) M 3/2 

1 1 
2nAv F ^2D(F + 1) |w| 3 / 2 
V2F 2 



(A.3a) 
(A.3b) 



7 rA I / FV / D( J F 1 + 1)C 



V2F 2 



7rA(F + 2) (VfTT+i) ^f|w| 3 / 2 v^D(FTi) 

(A.3c) 



HereC = (VV - ^1 + v A FTT V / M)(I w ~ + ( F + 
This result was also obtained in Ref. [5]. For d — 2 
the kernels are 



1 + 2F 2 + Sf 



ttD (4F 2 + 12F + 8) dis F \uj\ 
1 



8nD(F + l)dv F \u\ 



(A.4a) 
(A.4b) 



F 2 



2ttD(F + l)^ F (|c - F Ms | + (F + 1)| W |) (A 4c) 



2irD(F 2 + 3F + 2) dis F \uj\ 
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and for d = 3 we get 

IC a = t= F+2 — (A.5a) 

2tt 2 u f ^J\(2D{F +I)) 3 / 2 

K b = =i (A.5b) 

2n 2 v F J\^\{2D(F + l))3/2 



F(|(F+lM 3/2 -|"-*>»l 3/2 ) 
c V2n^ F (D(F + 1)) 3 / 2 (m s + w)((F + 2)w - f> s ) 
m.^o F ((F + l) 2 - VF+T) 
^ V2tt 2 D 3 / 2 (F + l)*(F + 2)v Fy /\u\ 

(A.5c) 

In the present paper, we only use the kernels with fi s — > 0. 
This is because in the limit // s C T the effect of a finite 
£i s is quadratic in fi s . 
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